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Abstract 

We study the roughness of randomly crumpled elastic sheets. Based on 

analytical and numerical calculations, we find that they are self affine with 

a roughness exponent equal to one. Such crumpling occurs e.g. when wet 

paper dries. We also discuss the case of correlated crumpling, which occurs 

in connection with flattening of randomly folded paper. 
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Take a sheet of paper. Wet it with water and let it dry. The sheet will crumple.Q A 
closer look at the dried surface does not reveal any preferred length scale in the crumpling. 
Rather, smaller crumples appear inside the larger ones, and so on to smaller and smaller 
scales, finally to reach the scale of the fibers making up the paper. Surfaces having this 
appearance are in fact abundant in nature. For example, fracture surfaces show roughness 
on all scales, as do surfaces that have been corroded |]-§]]. Recently, mathematical tools 
have been developed to describe such surfaces in terms of their scaling properties |5|||. We 
will in this letter analyse the roughness of dried paper surfaces in light of these scaling 
properties based on a theoretical model originating from the theory of linear elasticity of 
thin plates J7|||. We also note the connection between this problem and the problem of 
unfolding a randomly wrapped paper ||: Wrap a sheet of paper randomly into a tight ball. 
Then unfold it and try to flatten it as much as possible. It will stay rough by much the 
same mechanism that makes dried paper rough. However, the foldings of the paper have 
introduced long-range correlations between the local crumpling of the paper. These long- 
range correlations may change the roughness of the sheets. Experimental measurements 
of this roughness show, however, that it is within the experimental uncertainty equal to 
the roughness due to uncorrelated crumpling. We generalize the one-dimensional model of 
reference || to a two-dimensional elastic sheet. However, the roughness found in this model 
is larger than the one found in the random crumpling case and in the experiments. 

The scaling properties alluded to above, are more precisely described through the concept 
of self affinity. We choose for a given surface the (x, y) plane to be the mean plane, and the 
z axis to be the normal direction. The surface is self affine if it is (statistically) invariant 
under the rescaling transformation 

lr The baking of traditional Norwegian fiatbread or Indian papadam produces a much stronger 
crumpling, but of the same kind. 
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where A is a arbitrary scaling factor, and fi is a function of A. Group properties on the 
transformation (jl]) imply that fi is a homogeneous function of A characterized by an exponent 
H, fi oc X H . This exponent is the roughness index. The aim of this letter is to determine 
the roughness index of a randomly crumpled elastic sheet. 

We return to the problem of wetting and then drying paper. By pouring water on one 
side of a sheet of paper, it immediately curls up away from the wetted side. This is due 
to the paper expanding locally on the wetted side. If, however, the paper is wetted on 
both sides, a uniform expansion occurs, and the now thoroughly wet paper will stay flat. 
However, as it begins to dry, it roughens. There are two mechanisms behind this roughening: 
(1) internal torques which bend the surface locally due to inhomogeneous contractions in 
the perpendicular direction of the sheet, which in turn are due to inhomogeneities in the 
structure of the paper itself and/or differences in the drying drying rates from the two sides 
of the paper, and (2) buckling due to locally homogenous in-plane contractions of the paper. 
In deciding which of these two mechanisms is the dominant one, we note that buckling, 
which is a non-linear elastic problem, produce multiple stable bending modes, while the first 
process to a first approximation is a linear problem producing unique bending modes. Using 
whether a unique roughened shape of the paper results when the paper has dried, or there 
are several equivalent ones, we find that mostly the former is the case. We thus conclude 
that it is the internal torques which dominate. 

When wrapped paper is flattened, the resulting roughness is caused by internal torques 
and not buckling. 

In light of the above discussion, we base our modeling on the linearized equation for 
flexion of an elastic plate under a loading q(x, y) in the z direction: 

DAAC = q, (2) 
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where D is the flexural rigidity of the plate, (chosen to be unity, D = 1, in the following), 
and £(af) is the vertical deflection of the plate at x = (x, y) Equation (^|) is known as the 
Germain equation.^ However, the forces acting on the surface in the case of drying paper 
are not external. We are rather seeking a behavior where locally the surface at x would 
adopt a non-zero curvature A( = i](x) if it were isolated from its neighborhood. The 77 (x) is 
chosen here to be a random field. This noise term is assumed to represent the fluctuations 
in curvature due to the local inhomgeneous contractions of the paper. However, A( is not 
a force, AA( is, see equation (H). Therefore, the equation that balances the forces locally 
in the plane is 

AAC = At? . (3) 

With the assumption that r](x) is an uncorrelated random field, we have that (rj(xi)r](x2)) = 
c 2 S(xi — X2), where (. . .) signify ensemble average, and 5(x) is the Dirac distribution. Note, 
however, that modeling the unfolded and flattened paper, the stochastic field r\ has long 
range correlations build into it stemming from the folding process. 

Let us assume a system of size Lx L where L is large. The height fluctuations, ((C(^i) — 
((x 2 )) 2 ), may easily be calculated by using Green's method. We have that 

((x) = J dy G(x - y)A y r](y) = J dy A x G(x - y)r](y) , (4) 

where G(x) is the solution of AA£ = 5(x), which in Fourier space is G{k) = l/k A . Thus, we 
have that 



(C(fi)C(x 2 )) = j J dyidm A^G^fi - yi)A X2 G(x 2 - y 2) {r]{y 1)^2)) (5) 
= c 2 J A; 4 dkG(k)G(-k)e l{gl -* 2) * = c 2 J ^j^-^ . 

We may now calculate the height fluctuations: 



2 An interesting note on the history of this equation and Sophie Germain can be found in [10|. 
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<(CO?i) - C(^)) 2 ) = 2c 2 j o J — (l - e *l*i-<*l*«»W) = cm _^ 2 |2 ^ (6) 

where C is a constant. This last result shows that the roughness index is 

H = 1 - (7) 

One should also note that the variance of the height scales as 

<C 2 (*)> cc c 2 L 2 (8) 

with the system size. The square dependence on L is also a reflection of H = 1. 

A direct way to model numerically surfaces following equation (^), is to integrate this 
equation to get 

AC = v + <P , (9) 

where is a solution of the equation A0 = 0. We now assume that our system of size L x L 
has biperiodic boundary conditions. Then, we have that (A£) = if the field ( is to be 
continuous. Equation (|) then leads to the condition (<f>) = —(rj) on 0. The solution of the 
Laplace equation with biperiodic boundary conditions is a constant, <ft — constant = —(77). 
Thus, we have in Fourier space that 

C = ^(v-(v))- (10) 

This is straight forward to implement on the computer. 

In addition to generating the solution of equation (|3p through equation (TO), we have 
discretized (H) an d solved it through an iterative method. This method raises principal 
questions in that the noise term r\ is not a smooth function, but varies rapidly on all scales, 
including the one chosen as our lattice constant in the discretization. This is a typical 
problem faced when modeling strongly disordered systems. However, we note that a priori 
a continuum description of a strongly disordered system (where this is at all possible) is not 
more fundamental than a description based on a discrete system. 
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The discretization we have chosen is the following one: We use a triangular tiling of 
rhombical plaquettes as shown in figure |1] JET]] . The flexural energy associated with plaquette 
0394 (the numbering refers to the nodes of figure [I]) is £0394 = (Co — C3 + C9 — C4 + Vu) 2 , 
where ^34 corresponds to the noise 77 of equation (Q). Each elementary triangle is covered 
by three plaquettes, e.g. triangle 034 is the intersection of plaquettes 0394, 0345 and 0234. 
With this construction, the dependence of the elastic energy of the network on a given node 
in figure [l] is given by the sum 
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The total elastic energy of the network is the sum E = Eili e t /12, where i runs over all N 
nodes. In the limit of a vanishing lattice constant, 5, deo/d(o = leads to the equation 

5 4 AA( = -A V + —AAr ] + ... . (12) 

When 7] is a smooth function in the scale of 5, the second term on the right hand side of 
equation (|T2|) is small compared to the first term. However, if 1] is not smooth, it will be of 
the same order as A ~ l/<5 2 , and the second term (and those we have not explicitly written 
in (|T2|)) will be of the same size as the first. This is the case we deal with here. However, even 
though the implementation (|ll|) does not reproduce (|3]) on small scale when r] is not smooth, 
the large scale features will be the same as predicted by the continuum equations. Had our 
basic equations not been linear, however, such a conclusion could not have been drawn. We 
note, furthermore, that equation ( [11]) is an equally accurate but distinct modelization of the 
original problem, as the one based on the continuum theory, (^). 
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The noise term rjij in (|ll|) is associated with the bond ij, which in turn is associated 
with a plaquette. We choose the value of r/y to be either 1 or —1 with equal probability. 
We implement biperiodic boundary conditions. For a given distribution of noise terms rjij, 



we minimize the total energy E using a conjugate gradient algorithm as described in [JT2 
In figure |2|, we show a typical crumpled surface of size L x L = 256 x 256. 

In order to measure the roughness index we generated an ensemble consisting of 1000 
networks of size L = 10, 1000 samples of size L = 20, 800 samples of size L = 30, 600 
samples of size L = 40, 200 samples of size 60, 80 samples of size L = 80, 40 samples of 
size L = 100, 20 samples of size L = 150, and 10 samples of size L = 200. For each sample 
we measured the width defined as w 2 = (J2i(f)/N — (J2i G) 2 /-^ 2 ; which corresponds to 
and then averaged this quantity over all samples with the same size L. The result is shown 
in figure |[ We find that w oc L H with H = 1.0 ± 0.05, in agreement with the predicted 
H = 1. In figure |] we show the power spectrum P(k), which is the Fourier transform of the 
height-height correlation function, measured along one-dimensional profiles and averaged 
over 10 samples of size 256 x 256. k is the wave number. For a self-affine surface, one 
expects P(k) oc /c - ( 1+2i ^. We find a slope equal to —3.0 ±0.1 in figure f|, again consistent 
with a roughness index equal to one. 

Thus, we have succeded in demonstrating how a simple linear elastic theory of crumpled 
surfaces can generate non-trivial self-affine surfaces. It would be very interesting to see how 
these theoretical results match with experiments. We furthermore note that we may use 
these results to generate large self-affine surfaces by e.g. spotwise random heating of elastic 
plates. 

We now turn to the problem of randomly folded paper. The folding process induces in the 
paper a bending along the fold line somewhat similar to that observed when wet paper dries. 
Suppose now we fold the paper a first time. This produces a straight bend. However, the 
next fold — when the paper is unwrapped again — will change direction when it meets the 
first fold, unless this second folding is perpendicular to the first one. In addition, the bending 
of the second fold will change "sign" when it crosses the first fold. Thus, taking as starting 



point the one dimensional geometrical model introduced in reference 0, we implement the 
following two-dimensional model. We associate with each plaquette in a square lattice of 
size 2 N x 2^ a value 77. Assign to the line of plaquettes with coordinates (2W— 1), j), where 
j = 1, ... 2^, equal values for rj, either +1 or —1. Then, go to line (i, 2 Ar_1 ), where i — 1, 
2 N ~ 1 . Assign to the first node an revalue equal to either +1 or —1 with equal probability. 
This value is reproduced for all subsequent nodes along this line until the previous line is 
encountered at {i = 2 N ~ 1 , 2 7V ~ 1 ). Then, for all (i > 2 N ~ 1 , 2 iV ~ 1 ), we assign the opposite value 
for rj. We then go to the line (2 N ~ 2 ,j), choose an revalue equal to ±1 for node (2 N ~ 2 , 1), 
reproduce this value until we cross an older line, and then change sign for the subsequent 
rj values. Each time an older line is met, we change sign in the assigning of 77 values. This 
construction is repeated over and over until all nodes has a value for 77 associated with it. 
We show in figure |5| a pattern of bending lines produced this way for a 64 x 64 lattice. Black 
line means an associated f] = 1, and gray means that the associated 77 = —1. 

We then proceed to solve equation ([$]) with this distribution of 77 values. In practice, we 
used the solution given in equation fliTf). We show in figure |5] a 256 x 256 sheet produced this 
way. We note the similarity with figure ^ that was generated from a random distribution 
of r\ values. However, an analysis of the power spectrum of an ensemble of such lattices (20 
samples of size 256 x 256) gave the result shown in figure [7], where we plot on log-log scale 
P(k) vs. k in the same way as in figure |4j. There is a rich structure in this power spectrum 
which reflects the correlations induced by the folding process. Attempting to extract a power 
law from this structure, we show astraight line that follows the densest region of the plot. 
This curve has a slope of —4, thus indicating a roughness index equal to 1.5 ± 0.1 (from 
1 + 2H = 4). This is appreciably larger than what was found for the randomly crumpled 
paper, H = 1, and the experimental study of the randomly folded paper ||: Hi « 0.9. 

We have, thus, shown in this letter how a random internal bending process of an elastic 
plate leads to a self affine structure with roughness index equal to one. Such processes are 
met e.g. when wet paper dries. We have also studied the case of a correlated bending 
process, i. e. when paper is crumpled into a tight ball and then unfolded and flattened. We 
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find through a simple model for the correlations a roughness index which is higher than the 
one found in the random crumpling case. 
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Figure 1 

F. Tzschichholz, A. Hansen and S. Roux 

' 'Roughness of randomly crumpled elastic sheets ' ' 



FIG. 1. The numerical calculation is done using a triangular tiling of rhombical elastic pla- 
quettes as shown in this figure. 
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Figure 2 
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' 'Roughness of randomly crumpled elastic sheets ' ' 



FIG. 2. A randomly crumpled sheet of size 256 x 256. 
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Figure 3 

F. Tzschichholz, A. Hansen and S. Roux 

'Roughness of randomly crumpled elastic sheets ' ' 



FIG. 3. Average width w of sheets of size L x L. The slope of the broken line is 1.0. 
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Figure 4 

F. Tzschichholz, A. Hansen and S. Roux 

' 'Roughness of randomly crumpled elastic sheets ' ' 



FIG. 4. Power spectrum P(k) as a function of wave number k measured along one-dimensional 
profiles of surfaces that have been randomly crumpled. The curve is an average over all 256 profiles 
of each of 10 samples of size 256 x 256. The slope of the broken line is —3.0. 
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Figure 5 

F. Tzschichholz, A. Hansen and S. Roux 

' 'Roughness of randomly crumpled elastic sheets ' ' 



FIG. 5. A folding pattern generated by the algorithm described in the text. Black line means 
the bend is "upwards" while a gray line means that the bend is "downwards." The size here is 
64 x 64. 
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Figure 6 
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' 'Roughness of randomly crumpled elastic sheets ' ' 



FIG. 6. A sheet that has been been crumpled through the random folding algorithm. The 
size is 256 x 256. 
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Figure 7 
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' 'Roughness of randomly crumpled elastic sheets ' ' 



FIG. 7. Power spectrum P(k) as a function of wave number k for the randomly folded surfaces. 
We measured along one-dimensional profiles over 20 samples of size 256 x 256. The slope of the 
broken line is —4.0. 
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